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Band-inverted electron-hole bilayers support quantum spin Hall insulator and exciton condensate 
phases. We investigate such a bilayer in an external magnetic field. We show that the interlayer 
correlations lead to formation of a helical guantum Hall exciton condensate state. In contrast to the 
chiral edge states of the quantum Hall exciton condensate in electron-electron bilayers, existence of 
the counterpropagating edge modes results in formation of a ground state spin-texture not support¬ 
ing gapless single-particle excitations. This feature has deep consequences for the low energy behav¬ 
ior of the system. Namely, the charged edge excitations in a sufficiently narrow Hall bar are confined 
i.e. a charge on one of the edges always gives rise to an opposite charge on the other edge. Moreover, 
we show that magnetic field and gate voltages allow to control confinement-deconfinement transi¬ 
tion of charged edge excitations, which can be probed with nonlocal conductance. Confinement- 
deconfinement transitions are of great interest, not least because of their possible significance in 
shedding light on the confinement problem of quarks. 

PACS numbers: 


The role of the electron-electron interactions for the ex¬ 
perimentally accessible topological media is best appre¬ 
ciated in quantum Hall (QH) systems. The fractionally 
charged quasiparticles have been studied at the fractional 
filling factors mm, and the non-abelian excitations of 
more exotic QH states may eventually lead to a revolu¬ 
tion in quantum computing [3Hz!. However, Coulomb in¬ 
teractions play a crucial role also in the case of the integer 
filling factors [iiiHHia. Remarkably, interactions create a 
QH ferromagnetic ground state at = 1 even in the ab¬ 
sence of Zeeman energy. In such systems, the SU{2) spin 
rotation symmetry is spontaneously broken, resulting in 
the low-energy excitations being spin waves and charged 
topological spin textures, skyrmions [niainiiig. The 
presence of a small symmetry-breaking Zeeman field does 
not change the low-energy excitations qualitatively. 

In QH bilayer systems the role of spin is played by the 
layer index (pseudospin) [H [101 El [13] ■ In this case the 
SU(2) pseudospin rotation symmetry is explicitly bro¬ 
ken by the interactions, as they are larger within the 
layers than between the layers. The interactions favor 
the pseudospin orientations in (a;, j/)-plane, where the di¬ 
rection is chosen spontaneously (spontaneous 17(1) sym¬ 
metry breaking) so that the QH bilayers realize an easy- 
plane ferromagnet. Since the spontaneously chosen di¬ 
rection in the (x, 2/)-plane corresponds to a spontaneous 
interlayer phase coherence, this easy-plane ferromagnetic 
state is equivalent to an exciton condensate duns]. 

The QH ferromagnet and QH exciton condensate in 
electron-electron bilayers support a single chiral edge 
mode. However the two Landau levels may also support 
counterpropagating edge modes. The natural hosts of 
such kind of QH states are systems supporting quantum 
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Figure 1: Four different types of experimentally accessible QH 
pseudospin ferromagnetic states at iz = 1. The single layer 
realizations (left column) realize Heisenberg ferromagnets be¬ 
cause the interactions have SU{2) symmetry. The bilayer 
QH exciton condensates (right column) can be described as 
an easy-plane ferromagnet with a spontaneously broken U{1)- 
symmetry. We argue that the classification is additionally en¬ 
riched as the QH systems can support either chiral (top row) 
or helical (bottom row) edge excitations. 


spin Hall (QSH) effect [TIHTU] due to inverted electron- 
hole bandstructure. In these materials the magnetic field 
allows to tune through the Landau level crossing [20l [21] , 
where we expect to find a QH state with spontaneously- 
broken (pseudo)spin-rotation symmetry. Thus, we argue 
that there exist four different experimentally accessible 
pseudospin ferromagnetic states, determined by sponta¬ 
neously broken symmetry [SU(2) in single layer and U(I) 
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in bilayer systems] and the edge structure [chiral or heli¬ 
cal]. All these possibilities are illustrated in Fig. 

In this paper we concentrate on the helical QH exci- 
ton condensate state [broken U(l) symmetry and helical 
edge structure]. Remarkably, we find that in this sys¬ 
tem the charged edge excitations in a sufficiently nar¬ 
row Hall bar are confined: A charge on one of the edges 
is always connected to the opposite charge on the other 
edge through the bulk by a stripe of rotated pseudospins, 
and thus low-energy isolated charged excitations cannot 
be observed. The gapless single-particle excitations are 
prohibited since the electron-electron interactions lead to 
an edge reconstruction and opening of a single-particle 
gap [221 I2S]- However, unlike it happens in the exist¬ 
ing examples, the helical exciton condensate creates long- 
range correlations between edges. We show that a mag¬ 
netic field and gate voltages can be used to tune in and 
out of the exciton condensate phase. Thus this system 
provides a unique opportunity to study a confinement- 
deconfinement transition, similar to the one which is hy¬ 
pothesized to liberate the quarks from their color con¬ 
finement at extremely high temperatures or densities |24j . 
Finally, we show that the confined and deconfined phases 
can be distinguished using nonlocal conductance. 


I. HELICAL QUANTUM HALL EXCITON 
CONDENSATE PHASE 

We consider bilayer QSH systems, such as InAs/GaSb 
[I7HI9!, described by the BHZ Hamiltonian [niiiTiES]. 
The important property of these systems is that there is a 
crossing of electron and hole Landau levels as a function 
of magnetic field a,t B — Hcross (2011111 [25] [see Fig.j^a)], 
where the band inversion is removed. Near this crossing 
the single-particle Hamiltonian is |2S] 



Figure 2: Magnetic field and momentum dependence of the 
lowest Landau level energies, (a) There exists a robust cross¬ 
ing of the lowest Landau level energies as a function of mag¬ 
netic field at B — Rcross, because the energy of the electron¬ 
like Landau level with spin up (red line) increases and the 
energy of the hole Landau level with spin down (blue line) 
decreases as a function of B. We denote the energy separa¬ 
tion between these two levels as 2EGb- (b-c) Momentum (or 
equivalently position) dependencies of the lowest Landau level 
energies. The electron-like (hole-like) Landau level bends up¬ 
wards (downwards) in energy near the edge of the sample. 

(b) For B < Rcross the system supports helical edge modes 
protected by spin-resolved Chern numbers Cf = — C 4 , = 1. 

(c) For B > Bcross the edge is gapped according to the non¬ 
interacting theory (U-i- = = 0 ). 


to topological reasons Ea{y) > 0 reaches extremely large 
values (on the order of the energy separation between 
the bulk Landau levels) close to the edge |2S|. Therefore, 
for the magnetic fields B < Baoss, Eciy) goes through 
zero near the edge, yielding to the helical edge states [see 
Fig. [^b)]. On the other hand for B > Bcmss the edge is 
gapped in the non-interacting theory [see Fig. [^c)]. 

The electron-electron interactions E[j are described by 


Ho = - '4>l,i^k,i]EGikll), ( 1 ) 

k 

where EQ(y) = Eciklg) is the energy-momentum dis¬ 
persion of Landau levels and are the electron cre¬ 

ation operators for the lowest electron and hole Landau 
levels. Here we have fixed the total filling factor of the 
Landau levels vt = = 1; and utilized the fact 

that the momentum k in the Landau level wavefunctions 
is directly connected to the position y in the real space. 
Importantly, the spin and layer degrees of freedom are 
locked with each other, so that the pseudospin f (i) 
means simultaneously up (down) spin and upper (lower) 
layer [25]. The Fermi level is set to be at zero energy. 

In the bulk the energy Eg (y) = EGb is independent of 
the momentum {EGb < 0 for B < Baoss and EGb > 0 
for B > Hcross)- When approaching the edge the Landau 
level originating from the electron (hole) band always dis¬ 
perses upwards (downwards) in energy. The spatial vari¬ 
ation of EG{y) occurs within a characteristic length scale 
^0 ^ Ib, which depends on the details of the edge, but due 


Hl = ^'^Yl ^P'"'^<l)'^ka'^k',a''^k'+qa4k-qa, 

(T,(7' k,k' ,q 

( 2 ) 

where {k — k',q) is obtained by projecting the 
Coulomb interactions to the subspace generated by the 
wavefunctions of the lowest Landau levels [25]. Here, we 
assumed that the higher Landau levels are energetically 
separated from the lowest ones by an energy gap larger 
than the characteristic energy scale of the Coulomb inter¬ 
actions Vc" = e^/( 47 reeoZ_B). We find that this assumption 
can be satisfied with the material parameters correspond¬ 
ing to InAs/CaSb bilayers [26] . 

To find the ground state of the Hamiltonian H = 
Hq + Hi, we consider states where the local direction 
of the pseudospin h(r) (|h(r)| = 1) varies in space. 
Because the Hamiltonian is translationally invariant in 
the x-direction, we assume that h(r) is independent 
of X [22] ■ By further noticing that the y-dependence 
translates to a momentum dependence of the pseudospin 
hi{y) = hi{klg), we can express our ansatz for the ground 
state many-particle wavefunction as a Slater determinant 
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Figure 3: Phase diagram. For Eat < —2(Vo^^ — Vq) the sys¬ 
tem supports an uncorrelated helical QH phase. In this phase, 
there is no spontaneous interlayer phase coherence in the bulk 
{h^i, + hyi, = 0 ) and the system supports helical edge states 
(guaranteed by spin-resolved Chern numbers = il)- For 
li^Gbl < 2{V(^^ — V(f) the system supports a helical QH ex¬ 
citon condensate phase. In this phase, there is spontaneous 
interlayer phase coherence (h^i, hyi, 7 ^ 0 ), and the system 
supports exotic confined edge excitations (see below). For 
Eob > 2(Vo^^ ~ V(f ) the system is in a trivial QH phase, 
where the edge is fully gapped. Ecb can be controlled with 
gate voltages or magnetic field . Charged bulk excitations 
are gapped everywhere in the phase diagram (Flgap.s > 0 ). 
We have chosen Vq = 1/2. 


the three distinct phases of the system, which are sum¬ 
marized in Fig. For sufficiently large \EGb\, we see 
that = 1, meaning that only one layer is occupied. 
The phases h^b = 1 {uncorrelated helical QH phase) and 
hzb = — 1 {trivial QH phase) are topologically distinct 
from each other. For Ecb < —2(V(j^^ — Vf^) and hzb = 1 
the system supports helical edge modes (the spin-resolved 
Chern numbers are = ±1)- On the other hand, in 
the regime Eob > 2 (^ 0 ^^ — V^f) and hzb = — 1 , the edge 
is completely gapped (the spin-resolved Chern numbers 
are C^/i = 0). Between these two phases is the heli¬ 
cal QH exciton condensate phase, where \hzb\ < 1 and 
thus + hyf^ 7 ^ 0. In this phase the direction of the 
pseudospin projection onto the {x, ?/)-plane is determined 
spontaneously. Because the pseudospin in this system 
labels spin and layer index simultaneously, this phase 
has simultaneously spontaneous in-plane spin polariza¬ 
tion and spontaneous interlayer phase coherence. 

The bulk gap for single particle excitations ilgap.s can 
be calculated using Hartree-Fock linearization 


i^gap.s = 2^{EGb - 2V,^hzb)^ + MVo^^)^hl, + hi,). 


|lE'[h(A:Z^)]), where for each momentum k we create an 
electron with pseudospin pointing along h(fcZg). To com¬ 
pute the ground state, we need to minimize the energy 
functional for such kind of pseudospin texture [13]. For 
the energy functional we obtain |25| 

S = So - ^ j ^ V^^{k - k')K{kll)h,{k'll) 

k,k' ^ i—x,y 

+ U#(fc - k')hz{kll)hz{k'l%)\ -I- ^ EG{kll)hz{kl%). 

^ k 

(3) 

Here V^{q) = [-V;Hq,0) + V;Hq,q) + V^HhO)]/^ 
and Vp^{q) = Vp^{q,q)/^ are the interaction coeffi¬ 
cients, which characterize the anisotropy of the interac¬ 
tions within a layer and between the layers. 

We start by considering an infinite system. In this case, 
the pseudospin direction h(r) is spatially homogeneous. 
By minimizing the energy-functional (§ , we find that the 
pseudospin direction h;, is determined by the parameters 
EGb and [^- Here acts 

as an effective magnetic field preferring the pseudospin 
direction along —sgn{EGb)ez- On the other hand, the 
interactions prefer the pseudospin directions within the 
(x,y)-plane {V^'^ > l)f), and the energy cost to rotate 
the pseudospin so that it points along the ^-direction is 
proportional to ~ Thus, as a balance between 
these two competing effects, the direction of the pseu¬ 
dospin is tilted away from the {x, ?/)-plane, resulting in 


Additionally to the single particle excitations, the heli¬ 
cal QH exciton condensate supports collective excitations 
m- the neutral pseudospin waves (Goldstone mode) 
give rise to spin and counterflow charge superfluidity, 
and the lowest energy charged excitations are topolog¬ 
ical pseudospin textures, which carry fractional charge 
Here = (1 ± hzb)!^ are the pseudospin re¬ 
solved filling factors of the different Landau levels. The 
energy required to create these charged excitations is 
slightly lower than i?gap,s m- 

We point out that although hzb = ±1 are topologically 
distinct phases, the bulk gap for creating charged excita¬ 
tions never closes, when one tunes from one phase to the 
other by controlling EGb- This is possible because the 
pseudospin rotation symmetry, which protects the exis¬ 
tence of spin-resolved Chern numbers as topological num¬ 
bers, is spontaneously broken in the helical exciton con¬ 
densate phase. The interacting BHZ model for bilayers 
shows somewhat similar behavior also at zero magnetic 
field, where a trivial insulator phase can be connected to 
a quantum spin Hall insulator phase without closing of 
the bulk gap, because of an intermediate phase where the 
time-reversal symmetry is spontaneously broken |3()] . It 
is also experimentally known that the exciton condensate 
phase with = 1/2 can be smoothly connected to 

uncorrelated QH state with = 1 and = 0 in con¬ 
ventional QH bilayers [31j . Experimental investigations 
of InAs/GaSb bilayers in the QH regime [TH1[32| are con¬ 
sistent with this prediction, because no gap closing has 
been observed as a function of magnetic field. 
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Uncorrelated helical QH state (deconlincd phase) 


Helical QH exciton condensate (conhned phase) 
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Figure 4: Spin textures. (a),(b) Ground state spin textures for the uncorrelated helical QH phase and the helical QH exciton 
condensate phase, respectively. In both cases there exists a robust domain wall, where the polar angle of the pseudospin 
magnetization 6{y) rotates from tt to 9), along the y-direction. In the uncorrelated phase 06 = 0 but in the helical QH exciton 
condensate phase 06 yf 0, tt indicating spontaneous interlayer phase coherence in the bulk. In both cases the ground states are 
degenerate for all choices of the constant azimuthal angle (f) of the pseudospin magnetization. (We have chosen (f) = 0.) (c),(d) 
Charged excitations can be created by letting the azimuthal angle <()(a:) to rotate along the ^-direction. In a closed system 
(obtained by connecting the ends of the sample to form a narrow cylinder) (f>{x) must rotate integer multiples of 2ti. The 
energy to create such kind of excitation in the uncorrelated helical QH phase scales as 5E ~ Idw/L, and the elementary 
excitations, which carry charge ±e, can be created independently on the different edges. On the other hand, deep inside the 
helical QH exciton condensate phase 5E ~ V^'^WfL and the elementary charges are ±iz.|-e. Furthermore, the charged edge 
excitations are confined: a charge ±i/^e on one of the edges is always connected to the opposite charge on the other edge by 
a stripe of rotated pseudospins through the bulk, and thus isolated charges cannot be observed at low energies. The charge 
density obtained from Eq. is shown with red (positive charge density) and blue (negative charge density) colors. 


II. CONFINEMENT-DECONFINEMENT 
TRANSITION OF EDGE EXCITATIONS 


We now turn to the description of the ground 
state pseudospin texture hz{y) = cos[0o(2/)], hx{y) = 
sin[0o(2/)] cos(0) and hy{y) = sin[0o(y)] sin(())) at the edge. 
(The ground state will be degenerate with respect to 
the choice of (j).) As discussed above close to the edge 
Eciu) > 0 takes large values, because the edge states are 
topologically protected to exist at all energies between 
the lowest Landau levels and higher ones. Therefore close 
to the edge 0o{y) = n. On the other hand, in the bulk 
^ 0 ( 2 /) = Ob = a.TCCOs{hzb)- This means that there always 
exists a domain wall, where Oq rotates from tt to Ob- Al¬ 
though the existence of the domain wall is a robust topo¬ 


logical property of the system, the detailed shape of Oq {y) 
and the length scale Idw, where this rotation happens, 
depend on the details of the sample [21] ■ The ground 
states in the uncorrelated helical QH phase {9b = 0) and 
helical QH exciton condensate phase {9b h ^) il¬ 
lustrated in Fig. 1^ (a) and (b), respectively. It turns 
out that the existence of spontaneous interlayer phase 
coherence, which distinguishes the two different phases 
of matter, also has deep consequences on the nature of 
the low-energy excitations in this system. 

By using a Hartree-Fock linearization for the ground 
state, we find that the single particle excitations are 
gapped also close to the edge, and the magnitude of the 
energy gap is determined by the Coulomb energy scale 
. However, similarly as for the case of a coherent do- 
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Figure 5: Charged excitation localized at one edge of the sam¬ 
ple in the helical QH exciton condensate. The excitation can 
be visualized as a meron-antimeron pair, where the meron is 
localized inside the sample and anti-meron - outside. These 
excitations become the lowest energy charged excitations in 
sufficiently wide samples (otherwise similar cylinder geometry 
as in Fig. 1^, where it is energetically favorable to break the 
stripe connecting opposite charges at the different edges by 
a creation of a bulk meron. Depending on the profile Ecij/) 
near the edge, the lowest energy excitation is charge-neutral 
or having charge ±e. The charge is determined by the pseu¬ 
dospin orientation in the center of the meron. The energy of 
such excitation is of order of Coulomb energy 


energy 5E ^ V.^'^ldw/L in the uncorrelated phase and 
5E ~ V^^WjL deep in the helical exciton condensate 
phase. Here = \ ^qVp'^{q)q^l\ characterizes the 
cost of exchange energy caused by V^(x) [35]. 

In the uncorrelated phase these excitations have a 
charge ±e. By inspecting Fig. we notice that because 
the spin points along z-direction in the bulk, there is no 
rotation happening in the bulk. This means that we can 
choose separate fields </>i(a;) and 4’2{x) for the two edges, 
so that these excitations can be created independently 
on the different edges much as in graphene [35] . 

The situation is dramatically different in the helical 
QH exciton condensate phase. There, the elementary ex¬ 
citations in a closed system have a charge ±v-^e. More¬ 
over, as illustrated in Fig. a charge on one of 

the edges is always connected to the opposite charge on 
the other edge by a stripe of rotated bulk pseudospins. 
Breaking the bulk pseudospin configuration costs an en¬ 
ergy comparable to the Coulomb energy, and thus iso¬ 
lated charges cannot be observed at low energies. This 
means that this type of charged edge excitations in the 
helical QH exciton condensate phase are confined. 

It is illustrative to consider what happens to the exci¬ 
tations in the helical QH exciton condensate, when the 
width of the sample is increased. Namely, the excitation 
energy increases proportionally to the width of the sam¬ 
ple 5E ~ V.^^W/L and eventually for IT ^ L it becomes 
energetically favorable instead of having a large area of 
rotating spins between two edges to create a bulk meron 
(see Fig. [^. This resembles the physics of quarks, where 
the growing separation of a quark-antiquark pair eventu¬ 
ally results in the creation of a new quark-antiquark pair 
between them. 


main wall in QH ferromagnetic state in graphene |23j , the 
lowest energy edge excitations are not the single-particle 
ones. Namely, the ground state is degenerate with re¬ 
spect to the choice of (j>, and therefore in accordance with 
the Goldstone’s theorem the system supports low-energy 
excitations described by spatial variation of (j){r). Due to 
the general relationship between the electric and topo¬ 
logical charge densities in QH ferromagnets [Il[33|35| 


e 9 (/) clcos[0o(2/)] 
47r dx dy 


(4) 


these excitations also carry charge, which is localized at 
the edges of the sample. In this section we illustrate these 
excitations in a closed system obtained by connecting the 
ends of the sample to form a narrow cylinder with width 
W and circumference L [53] . 

We start by considering this kind of closed system, 
where W L (see Fig. [^. This geometry is topologi¬ 
cally equivalent to a Corbino ring, which has been exper¬ 
imentally realized for QH exciton condensates [34l l35] . 
Using Eqs. @ and @ we find that the lowest energy ex¬ 
citations correspond to rotation of (/>(cc) by 27r and carry 
a net charge within one of the edges [3^. They have an 


III. LUTTINGER LIQUID THEORY AND 
NONLOCAL TRANSPORT 


To predict experimentally measurable consequences of 
the charge confinement, we consider nonlocal transport 
in an open system. By considering the time-dependent 
field theory for the pseudospin for a reasonably narrow 
sample in the helical QH exciton condensate phase we 
arrive at an effective one-dimensional Hamiltonian [25] 


E[ = dx 


2WffT 


n(a:)^ 


Wpsb 


{d:^(j){x)f 


(5) 


where Psb = sir? 6b /is the pseudospin stiffness 

and F = e^/(167rZ^(Ug'^^ — Uq^)) describes the inter¬ 
layer capacitance per unit area, which is strongly en¬ 
hanced from the electrostatic value by the exchange inter¬ 
actions. The one-dimensional charge densities in the dif¬ 
ferent edges (labeled 1 and 2) ^ 1 , 2 ( 2 ;) = =F^(l+cos0f,)|^ 
are always opposite and determined by a single field 4>{x), 
highlighting the confinement of the charged edge excita¬ 
tions. The one-dimensional theory describes a Luttinger 
liquid, and the so-called Luttinger parameter K in the 
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convention used in Ref. [35], is given by 

... Ib (l + COS0b)2 

4sin0, ■ 

The Luttinger parameter in quantum Hall systems deter¬ 
mines the conductance for ideal contacts Gcf = Ke'^/h 
[351 [37|. Because the pseudospin waves are charge neu¬ 
tral in the bulk, conductance decreases with W as Gcf oc 
XjW . It is important to notice that in this system Gcf 
describes the conductance for a counterflow/drag geome¬ 
try, where opposite currents are flowing in the two edges. 
The helical QH exciton condensate phase does not sup¬ 
port net transport current as long as the voltages eV 
are small compared to hvir/W [35] . This automatically 
leads to a remarkable transport property that character¬ 
izes the helical QH exciton condensate phase. Namely, 
by considering a nonlocal transport geometry shown in 
Fig. Q (a), where a drive current is applied on one of 
the edges and a resulting drag current is measured on 
the opposite edge, we find that necessarily /drag = /drive 
at small voltages. This should be contrasted to the un¬ 
correlated helical phase, where the charged edge excita¬ 
tions are deconfined. In that case, one has two indepen¬ 
dent Luttinger liquid theories for the two edges [2l| , and 
therefore one expects only a weak drag current due to the 
Coulomb force acting between the charges. For W Ib^ 
we expect that this effect is negligible compared to drag 
current in the confined phase. 

Finally, to estimate the critical current /c, where the 
relation /drag = .^drive breaks down, we notice that the 
maximum voltage is determined by the gap eVmax ~ 
hvTT/W. By using reasonable estimates = Vj5^^/2, 
Ob = 7r/2, W = 201b, Ib = 10 nm, 

= 14 km/s [^ |31], we find Ic = GcfVmax ~ 0.1 nA. 

IV. SUMMARY AND DISCUSSION 

In summary, we have predicted the existence of a 
helical QH exciton condensate state in band-inverted 
electron-hole bilayers. We have shown that the counter- 
propagating edge modes give rise to a ground state pseu¬ 
dospin texture, where the polar angle of the pseudospin 
magnetization 9{y) rotates from the boundary value tt 
to the bulk value 6b along the direction perpendicular to 
the edge. Low-energy charged excitations can be created 
by letting the azimuthal angle of the pseudospin polariza¬ 
tion to rotate along the edge. Remarkably, in a suffi¬ 
ciently narrow Hall bar these charged edge excitations are 
confined in the presence of spontaneous interlayer phase- 
coherence {9b 7 ^ 0, tt): a charge on one of the edges always 
gives rise to the opposite charge on the other edge, and 
thus isolated charges cannot be observed at low energies. 
Moreover, we predict the possibility to control 9b with a 
magnetic field and gate voltages. This allows to study 
a confinement-deconfinement transition, which occurs si¬ 
multaneously with the bulk phase-transition between the 




Figure 6: Nonlocal transport to demonstrate the charge con¬ 
finement. (a) The transport geometry. A drive current is 
applied on one of the edges and the resulting drag current is 
measured on the opposite edge, (b) Due to the charge confine¬ 
ment in the helical QH exciton condensate phase, the drive 
current necessarily gives rises to an opposite drag current on 
the other edge. On the other hand, in the uncorrelated (de¬ 
confined) phase, there is only a weak drag current due to the 
Coulomb force acting between the charges. 


helical QH exciton condensate phase {9b ^ 0, tt) and the 
uncorrelated helical QH phase {9b = 0). 

The helical QH exciton condensate phase can be ex¬ 
perimentally probed using Josephson-like interlayer tun¬ 
neling and counterflow superfluidity [D Elia Ea [Ml IMl 
1381 [39| . Moreover, because the pseudospin in this sys¬ 
tem describes simultaneously both the spin and the layer 
degrees of freedom, the helical QH exciton condensate 
phase can also be probed using the spin superfluidity 
and the NMR techniques [Sj. Perhaps it is even possible 
to use local probe techniques to image the confinement- 
deconfinement transition and the confinement physics as 
illustrated in Figs. [^ and [^ Finally, we have shown that 
the charge confinement also gives rise to a remarkable 
new transport property. Namely, a drive current applied 
on one of the edges gives rise to exactly opposite drag 
current /drag = /drive at the other edge. 

Our results for the confinement of the edge excitations 
may also be applicable to the so-called canted antiferro¬ 
magnetic phase, which is predicted to appear in graphene 
[30] . Similarly to the helical QH exciton condensate state 
considered in this paper, the canted antiferromagnetic 
phase is characterized by a spontaneously broken U(l)- 
symmetry in the bulk and a single edge supports only 
gapped meron-antimeron excitations |40| . 
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The phenomena of confinement stemming from the 
particle physics models [24] has been studied also in con¬ 
densed matter systems I1IH151- However, we expect that 
the combination of the different techniques for probing 
the helical QH exciton condensate phase will provide a 
more intuitive understanding and new perspectives on 
the confinement physics. 

We also point out that InAs/GaSb bilayers is a promis¬ 
ing system for superconducting applications, and edge¬ 
mode superconductivity has already been experimentally 
demonstrated in the QSH regime H®] ■ In the presence of 


superconducting contacts, the helical QH exciton conden¬ 
sate may provide a new route for realizing exotic nonlo¬ 
cal Josephson effects and non-Abelian excitations, such 
as parafermions 1312]. 
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Supplementary material for ’’Confinement-deconfinement transition due to spontaneous symmetry breaking 

in quantum Hall bilayers” 

. EFFECTIVE MODEL FROM THE BHZ HAMILTONIAN 



B (T> 


Figure 7: Landau level fan for the InAs/GaSb bilayer. The figure shows the energies of the lowest 6 Landau levels on each side 
of the Fermi level (in terms of the minimum of the absolute value of energy as a function of magnetic field) for the Hamiltonian 
Q with characteristic values of the parameters for InAs/GaSb bilayers M = —6 meV, C = —1.4 meV, B = —78.3 eVA^, 
V — —18 eVA^ and A = 0.62 eVA [6]. We have used the bulk values for the g-factors = —14.7 and gh = —9.3. For 
B < Bcross the spin-resolved Chern numbers are C-f = —C_i = 1 (helical edge modes) whereas for B > Serosa the spin-resolved 
Chern numbers are Cf = Cj, = 0 (trivial insulator). 


We consider bilayer QSH systems, such as InAs/GaSb bilayers, described by the BHZ Hamiltonian lUH] 

AkyUQTy + {C -'D[(ka: - ^ f + k^] jeroTo 

( 6 ) 

where a and r are Pauli matrices in spin and electron-hole basis correspondingly, A4 describes the distance between 
the bottoms of electron and hole bands (in the inverted regime Ai < 0), C is the chemical potential, S,I? < 0 
(|;B| > I'D]) determine the effective masses for the electron and hole bands, 5e(/i) are the (^-factors for electron and 
hole bands, respectively, and the magnetic field has been written in the Landau gauge with Is = \Jh/eB being the 
magnetic length. This Hamiltonian describes an electron-hole bilayer, where the electron band in one of the layers is 
made out of s-orbitals and the hole band in the other layer is made out of p-orbitals, so that the tunneling between 
the layers (proportional to A) is odd in momentum. There exist two strategies for constructing this kind of bilayer 
system. The first possibility is to use two different semiconducting materials where the s-like electron band in one of 
the materials and p-like hole band in the other material are inverted such as InAs/GaSb bilayers [2H1!- The second 
strategy is to use in both layers the same semiconductor where the s-like electron band and the p-like hole band are 
close in energy, so that in a gated device one can reach a situation where the s-like electron band is active close to 
Fermi energy in one layer whereas the p-like hole band is active in the other. A promising approach to realize this 
possibility is to construct a bilayer in such a way that each layer individually supports the QSH effect [5]. Here we 
have neglected the spin-orbit coupling terms arising due to structural and bulk inversion asymmetry. In InAs/GaSb 
bilayers these terms are estimated to be very small [B]. The Landau level spectrum for InAs/GaSb bilayers is shown 
in Fig. In this material the parameters of the model can be tuned with the help of gate voltages and widths of the 
quantum wells. Here we have used characteristic values of the parameters [6] and the bulk values for the p-factors. 

The two lowest Landau level wavefunctions for this model are 

fe,„ = ^^(l^)(l(0) 0 0 0(l))L (7) 

where 4>{^) = Notice that spin and orbital degrees of freedom are locked, so that the pseudospin f 

(),) means simultaneously up (down) spin and upper (lower) layer. This locking is caused by the tunneling term 


H = 


{M - B[{k^ - + k^] ](JoTz+A[kx - ^)cr 

, gek-sB f , Qhg-BB 

H -,- ^^(ro -f tA H -^- cr^(ro - r^). 
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proportional to A, which is linear in momentum. Due to the existence of this term the electron-like Landau level 
with spin down (hole-like Landau level with spin up) couples to a hole-like (electron-like) higher Landau level, and 
as a result of this coupling these Landau levels are well separated in energy from the Landau levels given by Eq. Q . 
Within the subspace generated by these wave functions, the single-particle Hamiltonian is 

^0 = ( 8 ) 

k 

Here "0^ and are the creation and annihilation operators corresponding to the electronic states described by 

Eq. Q, and we have fixed the chemical potential so that the total density corresponds to one of these Landau levels 
being hlled and the other empty {vt = + = ^)- The Fermi level is set to be at zero energy. For an infinite system 

in y-direction, the energy Eq is independent on momentum and given by Ecb = M. — + S’‘+9h Due to the 

presence of an edge the Landau levels obtain an energy-momentum dispersion. According to Eq. Q the momentum 
is directly connected to the position y in real space, so that this energy-momentum dispersion can also be written as 
a position-dependent energy Ec^y) = EQ{klg). The Landau level originating from the electron (hole) band always 
disperses upwards (downwards) in energy, when approaching the edge. Moreover, because the edge states exist at all 
energies between the lowest Landau levels and the higher ones, close to the edge Ec{y) > 0 reaches extremely large 
values, which are on the same order as the energy separation between the bulk Landau levels. The spatial variation 
of Ec{y) occurs within a characteristic length scale Iq, which for clean sharp edge is given by Iq ^ Is- However, 
because the edge state velocity Vg = in the non-interacting theory is directly related to Iq, the edge roughness 

and disorder renormalize Iq upwards (vg downwards), and hence Iq can be considerably larger than Ib- The explicit 
form of the wave functions is given by Eq. Q only if Iq ^ Ib: but our main results are expected to remain valid even 
if this condition is not satisfied. 

Assuming that there is a band inversion at zero magnetic field {A4,B < 0), there is a crossing of the lowest Landau 
levels at magnetic field Baoss = Aif{^ — where the band inversion is removed. For B < Bcmss, we notice 

that Ec{y) = Eob < 0 in the bulk but Eo{y) > 0 close to the edge, yielding helical edge states. On the other hand 
for B > Hcross, Ea{y) > 0 everywhere, and therefore the edge is gapped according to the non-interacting theory. The 
change in the edge structures shows up in the spin resolved Chern numbers. Namely for B < B„oss the spin-resolved 
Chern numbers are = —C^ = \ (helical edge modes), whereas for B > Hcross = 0 (trivial insulator). 


ENERGY FUNCTIONAL 


Within the subspace generated by the lowest Landau level wave functions, the projected Hamiltonian can be written 
as iL = Hq + Hi, where the interactions are described by 




a,a' k,k' ,q 

Here, the projected Coulomb interactions can be written as 

rl/4 


vr'ik - C.,) = ^ 


y-{k-k' - q)l% 

Ib 


To simplify the expressions we assume that the quantum wells are very narrow so that 


(9) 


( 10 ) 




dTreeor ’ 


= U^^(r) = 


dTreeoVr^ -I- (P 


( 11 ) 


This assumption does not change the results qualitatively. However, one should keep in mind that quantitatively the 
energy scales associated with the interaction effects are overestimated, because the finite width of the quantum well 
would reduce the effective interaction strengths. 

To compute the energy for a pseudospin texture we follow closely the approach developed in Ref. [3 We assume 
that the components of the pseudospin hi{y) = hi{klg) are slowly varying with respect to Zs, so that we can express 
the many particle wave function as 


MHkil)]) 


n 


1 

72il - h,{kll)] 


[hM) - ^hy{kll)]^l^ 




| 0 ). 


( 12 ) 
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Here h‘f{kl%) = 1. The energy functional can be obtained by calculating 

E[hikl%)] = {^[hikl%)]\Ho + Hi\'f[hikl%)]). (13) 

Using the Wick’s theorem, we obtain 

E[hikll)] =Eo + J2EG{kll)h,{kll)-J2\v^ik-k')h,{kll)Kik'll)+ V^^{k-k')h,ikll)h,ik'll)], (14) 


k.k' 


where 


V^{q) = 


-VlHq,0) + VlHq,q) + Vl^{q,Q) 


Vc 


aV^l J-c 


i—x,y 


dy In 




3— (y/^s—9 *b)^/2 


J-oo 


and 


aV^l 

V^^iq)=^-Vl^iq,q) = 


_ 1tUU„ g ^ r dx 

ALb 


4'^ aV^L j _ 

Here the characteristic energy scale of the Coulomb interactions is Vc = /{A-Keei^lB)- 


(15) 


(16) 


MEAN FIELD SOLUTIONS IN THE BULK 


Before describing the pseudospin texture at the edge, let’s solve the ground state orientation of the pseudospin in 
the bulk. By assuming a homogeneous solution in the bulk, we obtain 


E[h] = Eo 


2 r 


2Trl^ 


EobK-V^^hi-V^^^hi + K) 


(17) 


where Eob is the bulk value of Eo{y) (in the inverted regime Ecb < 0 and in the non-inverted regime Ecb > 0), 


Vi =Y.^^iq) = Vc\ 


TT d 

2 Ib 


and 


U’' = E = Vc\Ey'f^‘iErk 


d 

_\/21b_ 


(18) 


(19) 


By minimizing the energy given Eq. with a constraint h? = 1, we obtain 


1 , 


Epb 




^Gb _ ^ 1 

Mm -1 < 2 


- 1 , 


0 *'0 

Epb ^ _-| 
2(yxv_vZ) < 


( 20 ) 


The other components satisfy = 1 — /i^f, and the energy is degenerate with respect to the rotations in the 

(x, j/)-plane. 

To estimate the single particle excitation gap in the bulk, we can construct a mean field Hamiltonian by Hartree-Fock 
linearization of the interaction terms. This way we obtain 


u , _ \^(JA ,/,t ') f ^Gb + m^b rria^b - imyb \ ( V’fc.t 

Vk,i) + i^Uyb -Ecb - m,bJ Ufc.l 


k 

vXY 


( 21 ) 


where m^b = —‘^Vihzb, rrixb = —‘^Vi^h^b and niyb = —2Vi^hyf,. The bulk gap for single particle excitations is thus 

Es.p,s = 2^{EGb - 2Vih,bi + + hi,). (22) 
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In addition to the crossing of the Landau levels as a function of magnetic field, several other conditions need to be 
satisfied in order to realize the helical exciton condensate phase: (i) The other Landau levels at the crossing point are 
separated in energy so that they are not excited, (ii) The densities close to the charge neutrality point can be obtained 
so that and can be controlled with magnetic field and gate voltages, (hi) The layer separation described by 
d/lB can be made sufficiently small to reach the exciton condensate phase, (iv) The temperature can be made small 
enough to reach the exciton condensate phase, (v) The disorder should not be too strong. 

(i) As demonstrated in Fig. ^the typical energy gap to higher Landau levels in InAs/GaSb bilayers is on the order 
of 10 meV, which is significantly larger than the gap opened by the exciton condensate order parameter. We also point 
out that although the large gap to higher Landau levels simplifies the theoretical analysis, it may not be necessary 
for the existence of exciton condensate state because the interaction effects actually tend to enhance this energy gap 
further. In fact, in GaAs bilayers the lowest Landau levels are separated from the higher ones by the Zeeman energy 
which is a rather small energy scale, but nevertheless the spin is fully polarized in the exciton condensate state [HHII]. 

(ii) In InAs/GaSb bilayers the densities close to the charge neutrality point can be experimentally reached in gated 
devices both in the absence o m na [i3] and in the presence [14] of the magnetic field. 

(hi) In GaAs quantum Hall bilayers it is known that the exciton condensate phase appears for d/lB ^ 1.8 [TSItl9) . 
Using the typical parameters of the experimental samples Emm] we find that this condition can be easily satisfied 
in InAs/GaSb bilayers. 

(iv) The temperature needs to be smaller than the energy gap opened by the exciton condensate order parameter. 
Moreover, since we are studying a two dimensional system the actual transition to the exciton condensate phase is a 
Berezinskii-Kosterlitz-Thouless transition [7|, and the transition temperature can be estimated to be on the order of 
Kelvin. The transition temperatures measured in GaAs bilayers are consistent with this type of estimate [18]. 

(v) The mean free path should be long compared to the coherence length. In quantum Hall exciton condensates 

the coherence length is on the order of Ib and therefore this condition is easily satisfied. However, the disorder 
plays also another role in this system. Namely, the vortices are charged and therefore they can be nucleated by a 
sufficiently strong disorder potential I201EI]. Although the detailed consideration of the disorder is not the subject 
of this paper, we point out that most of the phenomenology of the quantum Hall exciton condensate state survives 
at least approximately also in the presence of the disorder-nucleated vortices . 


DOMAIN WALL AT THE EDGE 


In order to describe the domain wall at the edge we first write hz{klg) = cos[6{klg)], hx{kP^) = sin[0(fc^^)] and 
hy{klg) = 0. [There is a degeneracy in {hx, hy)-plan.e, so that we can choose a specific direction arbitrarily.] The 
energy functional [Eq. can then be written as 

E[0{kll)] =J2EG{kll)cos[0{kll)]-J2\ypik-k')cos[0{kll)]cos[9{k'll)] + V^^{k-k')sm[0{kll)]sm[e{k'll)]\. 

k k,k' ^ ^ 

(23) 

By minimizing this with respect to 0{klg), we obtain 


tan{0{klg)) 


‘2Ek'yp^ik-k')sm[0{k'll)] 

2Efc' V^ik - k') cos[0{k'll)] - Ecikll) ■ 


(24) 


This can be solved numerically by iterations. 

Glose to the edge Ea{y) > 0 takes large values, because the edge states are topologically protected to exist at all 
energies between the lowest Landau levels and higher ones. Therefore close to the edge 0{y) = tt. On the other hand, 
in the bulk 0{y) = 9b = arccos(/iz(,). This means the existence of a domain wall, where 9 rotates from tt to 9b, is a 
robust topological property of the system, which is not sensitive to the details of the sample. 

The width of the domain wall is Idw max{/s, Iq}, where Iq is the length scale where Eciy) changes in the vicinity 
of the edge and Is is an intrinsic length scale, which is determined by the balance between the energy gain obtained 
by rotating 9 and the corresponding loss of the exchange energy. A rough estimate is obtained by expanding Eq. 0 
around 9b and noticing that the loss of exchange energy is approximately determined by (see Section []) 








(25) 


This leads to a simple estimate 




-Vl){l-2hl,)- EGbh.b 


I 


(26) 
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This length scale Is diverges at the phase transitions hz = 1, Eob = —and hz = —1, Eat = 2(Vf^^ — Vq^) 

indicating that 0(j/) approaches the bulk value in a power-law like fashion (outside the phase transitions it approaches 
the bulk value exponentially). Deep inside the helical quantum Hall exciton condensate phase Ig Ib- 


CHARGED EDGE EXCITATIONS 


We now consider excitations hz{r) = h^{y) + Shz{r), hx{r) = -^/l — hl{r) cos[(/)(r)] and hy{r) = -^/l — hl{r) sin[( 5 ()(r)], 
which are described by slow spatial variation of the pseudospin within the (x, j/)-plane [</'(r)] and small fluctuations 
[(5/iz(r)] of the z-component around the ground state configuration h^{y) = cos[do(2/)]- 
For this purpose, we first write Eq. (14) 


E[hikll)] = Eo[hz{kll)] -J2\ - k')K{kll)h,{k'll) 


k,k' ^ i—x,y 


(27) 


By expanding {i = x, y) 


h,{kll)h{k'll) = hl{kll) 


K{kll){k'll - kll)h[{kll) + -h{kll){k'll - kllfh'likll)... 


(28) 


we obtain 


E[h] - Eoihz] ^ 

f d2/| 

i—x,y i=x,y ^ 


- 2^/| J 

f 

+ X h^{r)\/^hx{r)\ 

^ i=x,y i—x,y '' 


- + 27TII J 

f d\- 

^ i—x,y i—x,y ^ 

(29) 


Here, we have transformed from momentum space k to real space y by using y = kl^, restored the possibility that 
h(r) may also vary in the cc-direction (using rotational symmetry of the system), assumed a boundary condition 


(fr V • {h,Whi) = 0 


(30) 


and is given by Eq. (25). 


Similar gradient expansion cannot be done for hz{r) because of the long-range Hartree interactions. However, deep 
inside the phases 6hz{r) mode is gapped by a large energy gap, and therefore we assume that as a first approximation 
it is enough to take the Hartree interactions into account only via the term. In order to calculate the excitation 
energy we define a functional 


E[h,X]=E[h]+^ J d\ A(r)(l-^/i,(r)2), 


(31) 


where we have now explicitly implemented the constraint /if = 1 with the help of Lagrange multiplier. The 
minimization of E [h, A] with respect to all variables leads to a saddle point equation (24) and allows us to consider 
fluctuations Shi{r) around the saddle point independently on each other. Eor general Ec^y) the saddle point equation 
can be solved only numerically. Here we proceed by assuming that Eciy) is slowly varying. This is always satisfied 
in the bulk and we will comment the influence of the edge corrections below. 

With these approximation the excitation energy (energy difference compared to the ground state) can be written 


as 


SE[Shz{r),(j){r)] = J + 




Eaiy) 


cos[6»o(2/)] 


(5/if(r) j. 
(32) 

A variational wave function describing a general pseudospin texture h(r) is no longer restricted to the lowest Landau 
level. Therefore, the many particle wave function needs to be projected to the lowest Landau level with an operator 
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7^, which can be implemented as explained in Ref. El This results in a charge density for such kind of excitations, 
which can be computed from the general relationship between the electric and topological charge densities [7] 


6p{r) = - —e^^h(r) • [9^h(r) x (9i.h(r)]. 

OTT 

For the low-energy excitations, the charge density is therefore (see also Section [ 


e d(j) d cos[9o{y)] 
Att dx dy 


(33) 


(34) 


From this expression we see that the charge density is localized in the vicinity of the edge, and the length scale in 
y-direction is determined by the domain wall width Idw Moreover, the low-energy charged excitations are always 
associated with spatial variation of (p along the edge, so that the charge density is proportional to dx4>- We will 
study the nature of these excitations separately in the two opposite limits Ecb ^ ~2(Vo^^ — (deep inside the 
uncorrelated helical quantum Hall state) and \Eci)\ 2(Vq^^ — Vff) (deep inside the helical quantum Hall exciton 
condensate phase). 


CHARGED EDGE EXCITATIONS IN A CLOSED SYSTEM 


We consider narrow quantum Hall systems with width W much smaller than the length L (similarly as shown in 
Fig. 4 in the main text), and we assume that the system is closed by connecting the ends of the sample. Based on 
Eqs. ( [^ and p4| we argue that the lowest energy charged excitations are obtained by letting azimuthal angle 
rotate along the x-direction. In a closed system the angle (j){x) needs to rotate integer multiple of 2tt. The lowest 
energy excitations in a closed system, which carry a net charge within one of the edges correspond to rotation of (l>{x) 
by 27r, and they have an energy SE ~ V^^ldw/L deep inside the uncorrelated phase and SE ~ V^^W/L deep inside 
the helical exciton condensate phase. 

In the uncorrelated phase we obtain from Eq. (34) that these excitations have a charge ±e. Moreover, they are 


deconfined: One can create these excitations independently on the different edges. 

In the helical quantum Hall exciton condensate phase, we obtain from Eq. (34), that the elementary excitations in 


a closed system have a charge is where = {1 + cosdt,)/2. Moreover, a charge Ev^e on one of the edges is 

always connected to the opposite charge on the other edge by a stripe of rotated pseuspins through the bulk. Creating 
this stripe costs an energy 5E ~ which is proportional to the distance W between the charges. However, 

breaking the stripe somewhere in the bulk would cost a much larger energy comparable to the Coulomb energy. Thus 
isolated charges cannot be observed at low energies. This means that in the helical quantum Hall exciton condensate 
phase the charged edge excitations are confined. 


LUTTINGER LIQUID THEORY 


To develop time-dependent theory for these excitations, we use the known result that the Euclidian action in the 
adiabatic approximation can be written as [7] 

S'^[h, A] = J ^''”1 y y A.(h) ■ 

where Vh x A(h) = h. Thus by assuming that Ec{y) is slowly varying function, we can expand the action around 
the saddle point and obtain 



SS^[P,6h,] 


1 

27rZ| 





sin^ do 



V(f)sin2 do(y) 


Eciv) 

2 


cos do(2/) 



(36) 


A. Confined phase 


Deep inside the confined phase, all the properties of the system will be determined by the bulk (see below). 
Therefore, we neglect the effects arising in the vicinity of the edge, and use Eq. (20) to rewrite the action as 


5S^[P,5K] = 


1 

2ttP^ 


dr / 


UhE 

2 


+ U/ Ij, sin^ do + (v - yo )SK 


(37) 
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We integrate out the massive fluctuations 5h^. This way we obtain 

By going to real time, we identify the Langrangian density as 

Ps{y), 


2e2 




( 38 ) 


(39) 


Here Ps{y) = sin^[0o(?/)]/7r is the pseudospin stiffness and T = e^/(167rZ^(V(j^^ — Vq)) describes the interlayer 

capacitance per unit area, which is strongly enhanced from the electrostatic value by the exchange interactions. (As 
an important check we notice that if one neglects the exchange contributions in , the expression for T 

simplifies to the familiar electrostatic formula for the interlayer capacitance.) 

The corresponding equation of motion is 


■■ 

V [ps{jj)v<i)\ = 0. 

This equation is translationary invariant in x and t thus we can express the solution as 

<^ = <^o + ^e-*[“"W‘-'="l/n.fc(y), 

n.k 


(40) 


(41) 


where (po is an arbitrary constant and fn,kiy) describe the different transverse modes with energy-momentum disper¬ 
sions En{k) = HuJnik). By substituting the trial solution (41) to Eq. (40), we And that the eigenfunctions fn,k{y) and 
eigenfrequencies uJn{k) satisfy an eigenvalue equation 


ti^Tujniky 


fn,k{y) + k‘^Ps{y)fn,k{y) - Ps{y)fn,k{y) - p's{y)f'n,k{y) = o. 


(42) 


which can be solved numerically. 

Numerical results for the helical quantum Hall exciton condensate phase show that the lowest energy mode fo,k{y) 
can be approximated as a constant throughout the sample. Moreover, the lowest energy mode at reasonably small 
momentum is separated from higher modes by an energy AE ^ hvTrjW. We restrict the analysis to sufficiently 
low-energies that these higher modes are not excited. Then we assume that cj) is constant in y-direction and the 
dispersion relation is 


uj{k) = 


e'^Psb 

h?r 


k = vk, 


where p^i, is the bulk value of Paiy)- Tbe Hamiltonian of the system then becomes 

Wpab 


H = dx 


2Wh^r 


n(x)^ 


{da:4>{x)Y 


(43) 


(44) 


where n(a;) is the momentum conjugate to (/>(a;). The one-dimensional charge densities in the different edges (labelled 
1 and 2) are 


PiAx) = T^(l + cos6»{,) —, 
47r ox 


(45) 


where 0b is the bulk value of 0{y). The confinement of the edge excitation clearly shows up here as the the charge 
densities on the different edges are always opposite and determined by a single field (p{x). This Hamiltonian describes 
a Luttinger liquid. In the standard convention [27] the Luttinger liquid theory is written as 


H = — dx 
2Tr I 


vKtt'^-, V 7,5 


ed(f) 

p{x) = -—. 

TT OX 


(46) 


To rewrite our Hamiltonian in this convention, we apply the transformation preserving canonical commutation rela¬ 
tions 


Pix) = 


p{x), n(a;) = ^ 


1 -I- cos 6b 


4 


(47) 
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and notice that the charge density at each edge is ^ 1 , 2 ( 2 ;) = ^p{x). 

This way, we can identify the Luttinger liquid parameters K and v as 

le^Psb _ e /l + coseb^ Ib (1 + cosdf,)^ 

’’“Vw “wy—lyr-<“*> 

Here, v is just the pseudospin wave velocity obtained earlier from the dispersion uj{k). The Luttinger parameter K 
determines the conductance Gcf = Ke^ jh. However, it is important to notice that Gcf describes the conductance for 
a counterflow/drag geometry, where opposite currents are flowing in the two edges. The helical quantum Hall exciton 
condensate phase does not support net transport current as long as the voltages eV are small compared to hvir/W. 


B. Deconfined phase 


In the uncorrelated helical quantum Hall phase the transverse modes are localized within the length scale Idw from 
the edges. Therefore, in contrast to helical quantum Hall exciton condensate phase we can define independent fields 
^ 1 ( 2 ) (cc) and charge densities y 0 i( 2 )(a:) = Thus charges can be created independently on the different edges, 

highlighting that in the uncorrelated helical phase the charged edge excitations are deconfined. 

We can now repeat the calculation done above for the confined phase. This way, assuming that Eciy) is slowly 
varying and that the transverse modes are described by a constant within a distance Idw in the vicinity of the edge 
and 0 elsewhere, we arrive to Luttinger liquid Hamiltonian 


^ = E 


dx 


2ldwh^^a 


-Xli{xf- 


^dwPs,Q 




(49) 


where ps 


fo'“dy V2^^sin^[6io(2/)]/7r and 


r — 

-L n/fj — 


dy eV{167rZ|[(t/o-^^ - Vjf) sin^ 6»o(y) - 


cos 00 (y)]}- This way the Luttinger liquid parameters are identified as 


Ps,av 


IL L n 


K = 


\/d av P s^o 


(50) 


However, in contrast to the helical quantum Hall exciton condensate phase the Luttinger liquid parameters depend 
strongly on the detailed shape of the Ea{y) function, and if it is not slowly varying in the vicinity of the edge, there 
can be large corrections to the expressions above. Nevertheless, the structure of the Luttinger liquid theory (491 is 
very general, and for physically reasonable parameters iL ^ 1. Thus we expect this phase to support helical edge 
state transport with conductance G ~ jh. Here G describes the conductance for a transport geometry, where a net 
transport current is flowing along one of the edges. 


. SPIN TEXTURE-CHARGE DENSITY RELATION IN QUANTUM HALL SYSTEMS AND 
NUMERICAL CALCULATION OF THE ELECTRIC CHARGE DENSITY FOR THE SPIN TEXTURES 


It is known that both topological and non-topological contributions to the electric charge exist at the vortices in 
various systems [281 - 135] . It is possible to show that the mechanisms considered in these systems do not give rise to a 
non-topological contribution to the electric charge in quantum Hall exciton condensates. However, instead of going 
through the mechanisms one-by-one, we present general arguments for the topological spin texture-charge density 
relation and discuss the conditions for its breakdown. We complement the analytical arguments with a numerical 
calculation of the electric charge density for the pseudospin textures considered in the paper. 

Originally the spin texture-charge density relation for quantum Hall systems was argued from the Chern-Simons 
relation between the density and the statistical magnetic field [36] . Namely, the electrons feel the spin texture via the 
additional Berry’s phase and therefore the orbital degrees of freedom are influenced in the same way as if additional 
magnetic flux density was inserted into the system. The extra magnetic flux is associated with an extra charge yeilding 
to Eq. (331 |36l [37]. This argument is valid only if the spin textures are smooth on the scale oi Ib- There exist also 


an alternative argument for the special case where the filling factors are = 1/2 and one neglects the higher 

Landau levels, so that there is a particle-hole symmetry in the system. In that case there is a general argument that 
a localized charge bound to a defect must be a multiple of ±e/2 [35]. Therefore, the possible charges ±e/2 appearing 
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in the special case = 1/2 can be explained this way (SHI- If the higher Landau levels are excited the particle 

hole symmetry is no longer exact. Therefore, this argument again relies on smoothness of the spin textures. Finally, 
the Eq. (331 can be obtained with an explicit microscopic calculation [7] and it can be demonstrated for specific spin 
textures with the help of explicit construction of the many-particle wave functions |39j . These approaches also rely 
on the assumption that the spin textures are smooth on the scale oi Ib- 

Based on these arguments it is clear that the relation between the pseudospin texture and the charge density 
[Eq. (33l] is valid if the pseudospin textures are smooth on the scale of Ib- The pseudospin textures considered in the 
main text [see Fig. 4 in the main text] satisfy this requirement. In fact we can explicitly construct the many-particle 
wave functions to illustrate how the charges appear in these textures. Namely, consider a hnite system with length 
L and width W illustrated in Fig. 4 in the main text. Because the momentum k in the Landau level wave functions 
is directly connected to the position y in the real space y = kl^, the possible values of k are kn = n27r/L, where 
n = 0,1, ...,N and N = WL/{2ttI^). The ground state pseudospin texture is given by h^{y) = hz{kl^) = cos[0o(fc^B)], 
hx{y) = hx{kl%) = sin[0o(A:/^)], hy{y) = 0 so that the many particle wave-function for the ground state |4 ')gs can be 
written as 


N 


i^)gs= n 


1 


„=o \/2{l-cos[6»o(A:„l|)]} 


sin[6lo(fc„Z|)]V]I^_^ + [1 - cos[6»o(/c„;|)]]i/|;^_^ ^0). 


(51) 


The elementary excitation is described by a pseudospin texture hz{y) = hz{kl^) = cos[9o{klg)], hx{x,y) = 
hx{x,klg) = sin[0o(fcZ^)] cos[(/)(a:)], hy{x,y) = hy{x,klg) = sin[0o(fc^B)] sin[( 5 ')(a::)], where (j){x) = 27rx/L. The winding 
of 4>{x) can be removed by introducing a momentum shift Ak = 27r/L between the electron and hole Landau level 
wave functions. Because the pseudospin-texture is slowly varying and the pseudopin close to the edge points down 
the many-particle wave function for the excited state can be written as 

^-1 1 r ^ ^ I 

I^)es = n “7or/ - “ cos[0o(fc„Z|)]][jO). (52) 

n=0 COS[(7o(,'SntBjjl L J 


From these expression one finds that the charge density is 

p r stt 

B V 


e d(j) dcos [6o{y)] 


47r dx 


dy 


(53) 


in agreerement with Eq. (341. From this expression it straightforwardly follows that the charges appearing at each 


edge for the pseudospin texture shown in Fig. 4 in the main text are where is the bulk filling factor for 

pseudospin up. It is also clear that in a closed system the charges of the possible edge excitations and the charges of 
the vortices must be related. As shown in Fig. 5 in the main text it is possible to end the charged edge excitation 
into a bulk vortex, and the total charge in the system must be an integer multiple of e. This requirement is satisfied 
because the possible charges of the vortices are and + vi = 1. 

In all arguments so far the whole pseudospin texture was considered to be smooth on the scale of Ib- However, it 
turns out that it is possible to relax this assumption. Namely, if one considers the type of pseudospin textures shown 
in Fig. 4 in the main text it is actually enough that the pseudospin texture is smooth on the scale of Ib in the bulk 
region between the localized charges. Close to the edge it may vary arbitrarily sharply. To understand this consider 
first the type of smooth pseudospin texture shown in Fig. 4, such that charge appears at one edge and —u^^e at the 
other. We can now deform one end of the pseudospin texture in such a way that the rest of the pseudospin texture 
remains smooth. Then, there is still charge z/^e localized on one edge and the bulk is charge neutral, which means that 
a charge —n^e necessarily remains at the other edge although it is no longer smooth on the scale of Ib- The values of 
the polarization charges in this system are therefore fully topological. No local perturbation in the pseudospin texture 
can change the charges appearing at the edges. Only a global perturbation which connects the fractionally charged 
excitations may give rise to redistribution of the charges and modify their values from ±z^^e. 

We have numerically verified the statements based on the analytical arguments given above. For this purpose we 
have inserted an order parameter term 


Ha = -Ao 


hz{x,y)azTo + hx{x,y)axTx + hy{x,y)axTy 


(54) 


to the BHZ Hamiltonian [Eq. ([^]. The charge density associated with the excitation 5p{y) can be numerically 
computed by calcuting the difference between the charge densities for pseudospin textures hz(x,y) = cos[0(?/)], 
hx{x, y) = sin[0(j/)] cos[(/(a;)], hy{x, y) = sin[0(?/)] sin[(/(a;)], where (j){x) = 2 'kx/I^ (excited state) and (j){x) = 0 (ground 
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(a) (b) 



Figure 8: Charge of the edge excitation in a closed system. Charge concentrated close to the edge of the sample as a function 
of (a) Ig and (b) l^. Here Ig determines the length scale where the angle 6{y) is varied from tt to 6b in the vicinity of the edge 
and I 4 , describes length scale where the angle (j>{x) rotates by 2tt. The charge q is computed by integrating the charge density 
over a single period along the x-direction and over the interval of width 51 b from the edge of the sample in y-direction. We 
have chosen Og = 7r/2 so that the expected charge shown by a dotted line is q = = e/2. As expected based on the analytical 

arguments the value of q is independent on Ig (as long as the width of the integration interval in y-direction is chosen to be 
larger than the width of the domain wall where the charge is concentrated in the vicinity of the edge), and deviations from the 
value e/2 as a function of start to occur only when l^i, < Hb- In the numerics we have chosen the order parameter strength 
Aq = 1 meV, the width of the sample W = 201b, and the space is discretized using a lattice constant 0.251b- For (a) we fix 
l(f, — 51b and for (b) Ig = 31b. The parameters of the BHZ Hamiltonian are same as in Fig. 


state). Here we have introduced a length scale where the phase <p{x) rotates over 27r. The charge of an elementary 
excitation q is obtained by integrating the charge density over a period l^f, in x-direction and over the width of the 
domain wall in y-direction. If we choose ^ Ib the charge of the excitation according to our analytic arguments 
should he q = ±z/^e. However, l^f, allows us also to control how fast the pseudospin direction changes in such a way 
that it influences the system everywhere in the bulk. Therefore, if we choose ^0 < it is possible to excite the higher 
Landau levels everywhere in the bulk so that topological protection is destroyed, and we can study how this affects 
the charge of the excitation. Additionally we can illustrate the topological nature of the charge of the pseudospin 
texture by demonstrating that deformations of the pseudospin texture appearing only close to the edge do not affect 
the value of q. For this purpose we define 

{ n + {9b- Tr)y/l 0 , 0 < y < Ig 

9b, le<y<W-lg (55) 

9b + {tt - 9b){y - W + le)/lg, W - Ig < y <W 

in such a way that Ig allows to control the length scale where the angle 9{y) is varied from tt to 9b in the vicinity of 
the edge. Based on our topological arguments the charge q should be independent of Ig. (We assume W Ig.) 

The numerical results are shown in Figs. As can be seen from these figures for smooth pseudospin textures the 
numerical results are in agreement with the analytical expectations. Moreover, q does not depend on Ig [Fig. ia)] 
in agreement with the topological argument. On the other hand, if we deform the pseudospin texture in such a way 
that it rotates fast in the bulk by choosing Hb it is possible to excite the higher Landau levels everywhere in the 
bulk so that topological protection is destroyed. This results in redistribution of the charges giving rise to deviations 
of q from the value [Fig. j^b)]. 
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